A chaotic 'turnstile' for atoms in periodic potentials 
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Using a new type of chaotic ratchet generated by pulsed standing waves of light, we propose a 
mechanism which would allow packets of atoms travelling through a pulsed optical lattice in one 
direction to pass almost undisturbed, while strongly heating atoms drifting through in the opposite 
direction. An analytical formula for the diffusive energy growth is derived and shown to give good 
agreement with numerical calculations. 
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There is much ongoing interest in ratchets- in other 
words the use of periodic (but spatially and/or tempo- 
rally asymmetric) systems which can be used to generate 
a current when there is no net force: see |jj for a com- 
prehensive review. Much of the work has focussed on 
Brownian ratchets Q which play an important role in 
biophysical systems such as molecular motors. Corre- 
sponding quantum ratchets have also been investigated 



There has been comparatively little work on Hamilto- 
nian ratchets, which are 'clean' ratchets without extrin- 
sic noise or dissipation. Mixed phase-space mechanisms 
(involving tori/stable islands) for ratchet transport have 
been investigated || [| . In ||, || general properties of 
ratchets and a mechanism involving desymmetrization of 
ballistic flights was proposed. In jfj a sum rule connect- 
ing the current in regular islands with that in adjoining 
chaotic manifolds was derived. 

Until recently it was not thought that a fully chaotic 
system could admit a ratchet effect. But in M a fully 
chaotic ratchet mechanism was proposed for atoms in 
a double-well periodic lattice, pulsed with unequal peri- 
ods. The ratchet depends on the generic properties of 
the chaotic diffusion so the initial state requires no spe- 
cific preparation: in all cases it was found that in general 
there is a classical timescale t r , during which particles 
moving with negative momenta (relative to the initial 
value) absorb energy at a different rate from those mov- 
ing in the opposite direction. In the corresponding quan- 
tum ratchet, for a significant effect to develop, it was 
found in || that it is simply necessary that the quantum 
break-time t* be of similar magnitude to t r . 

Atoms in pulsed optical lattices have become a 
paradigm in the study of quantum chaos; experiments on 
sodium and cesium atoms ^] have provided a convincing 
demonstration of dynamical localization (DL) |HJ [TT| , 
the quantum suppression of chaotic diffusion. The dy- 
namics for the usual experiment is given by the kicked- 
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particle Hamiltonian : H = ^ — K cos x^2 n S(t — nT) 
where K is the kick strength. The classical dynamics is 
obtained by iterating the well-known 'Standard Map' : 
Xi+i =Xi+ piT; p i+1 =Pi + Ksinx i+1 . 

We can take T = 1, without loss of generality, in the 
Standard Map, but for the ratchet, we use a repeating 
cycle of unequal kicks. The ratchet Hamiltonian is given 



by H = \ + V(x) J2 n j 5{t - nTj) where for the n-th 
cycle, the period of the j — th kick is Tj. For example, 
in H the kick spacings cycle between T\ = 1 + b, T2 = 
1,T 3 = 1 - b where \b\ < L 

In H an asymmetric potential of the form V(x) = 
K{sinx + asin(2ir + $)} was investigated, but in |7j] the 
spatial symmetry was broken by a 'rocking' linear term: 
V(x) = -(A'cosrr + Axsj) where Sj = . In @ the 

kick spacing alternated between Ti = 1 + 6 and T = 1 — b. 

The rocking ratchet was analysed in the regime 
where regular tori are still present ; it was found that al- 
though it no longer has the 2tt periodicity in p of the Stan- 
dard Map, there is a long range periodicity in p which is 
of order 2n/b. Also, since pairs of corresponding tori are 
located asymmetrically about p = 0, the resulting con- 
finement of classical trajectories can yield transport. 

Here we investigate this system in the chaotic regime. 
One aim of this work is to demonstrate that the physics 
is similar to that found in [|| and hence that the new 
chaotic ratchet mechanism is quite generic: both ratchets 
depend on the ratchet timescale t r ~ m t ne classical 
case, where D is the diffusion rate, and the break-time 
t* ~ D/ti 2 in the quantum case. The ratchet currents 
were found in ^ to originate predominantly from corre- 
lations of the form < V (xi)V (xi+2) > ■ We show below 
that the resulting analytical formulae accurately predict 
phenomena such as current reversals, without any de- 
tailed consideration of the structure of phase-space. 

Further, while [H considered only the ratchet current 
which arises for a system with zero initial current, ie 
< p{t = 0) >= 0, here we present results for the case 
where there is an initial drift < p(t = 0) >= pq. We 
show that this can form the basis of a device to manip- 
ulate traffic of cold atoms moving along some channel 
in a trap. This is not so far-fetched, since the trapping 
of atoms in devices such as atom chips is now a real- 
ity. The proposed device represents a sort of 'turnstile' 
which would selectively heat only atoms moving in one 
direction, while imparting little or no energy to atoms 
moving in the opposite direction, since we can tune the 
parameters of the pulsed lattice so as to generate corre- 
lations which almost cancel the energy diffusion in one 
direction, but enhance it in the other. Since the rock- 
ing ratchet produces much simpler analytical expressions 
than the double well ratchet, it provides a better didactic 
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example of the underlying physics. 

For the Standard Map , at the lowest level of approx- 
imation, the momenta at consecutive kicks are uncorre- 
cted and evolve in time as a 'random- walk'. The average 
momentum of a large ensemble of particles is unchanged. 
The average energy grows linearly: if the momenta are 
uncorrelated the average kinetic energy grows by K 2 /4 at 
each consecutive kick. In the absence of phase-space bar- 
riers the energy is unbounded and this diffusive increase 
in energy continues indefinitely. It is characterized by a 
diffusion rate D , ie < p 2 >= D t so for uncorrelated 
momenta Do = K 2 /2. However, this results neglects 
correlations between sequences of consecutive kicks; if 
included, they result in well-known corrections to the 
diffusion constant in the form of Bcsscl functions hence 
D = ^ {\-2{J l {K)) 2 -2J 2 {K)...) These cor- 

rections have even been measured experimentally with 
cold cesium atoms in pulsed optical lattices [ fL4| . 

Here, the C(2) = —K 2 J 2 (K) term is of particu- 
lar interest. It corresponds to a two- kick correlation, 
given by evaluating the phase-space average of 2 < 
V'(xi-i)V'(x i+ i) > hence: 

C(2) = 2K 2 < sinxi_i sina; i+ i > . (0.1) 

It is easy to show, by direct substitution from the map 
that one obtains separate integrals over the x and p . For 
chaotic diffusion it is usual to neglect odd integrals, for 
example < sinpcosp >~ 0, while even integrals yield a 
positive average, eg < sin 2 p >=< cos 2 p >= 1/2, since 
after a few kicks the ensemble has a substantial spread 
in p. In effect the only terms which contribute to C(2) in 
the standard map are even in p and do not differentiate 
between ±p. The integral over x yields J 2 (K) , hence the 
formofC(2). 

We consider now the map of (t]], with a rocking linear 
potential, which consists of a repeating cycle of two kicks. 
Around the i — th kick we have: 

Xi = Xi-x + + b) 

Pi = Pi-i + K sin Xi + A 
x t+ i = Xi +pi(l-b) 
Pi+i =pi + K sinx i+ i - A 

Following the usual procedure, we now obtain a modified 
form for the 2-kick average, but find that the momentum 
averages now include terms of the form < sin 2 p. sin 2pb > 
and < cos 2 p. sin2p6 >. At small t, pb is small while 
sin 2 p oscillates rapidly, hence we can approximate these 
by 1/2 < s'm2pb p av {t)b and 1/2 < cos2pb >~ 1/2 
where p av is the average momentum relative to the intial 
momentum. We note that such terms are the origin of 
the ratchet effect proposed in ||. These terms are odd 
in p and hence- unlike the standard map case- depend on 
whether we average our momenta from — > oo or from 
— > — oo. As in B below we consider the case of par- 
ticles with positive momenta separately from those with 



negative momenta, since the chaotic ratchet depends on 
differential energy diffusion rates for particles moving in 
opposite directions. 

We consider first the case where p(t = 0) = (or 
we take a cloud of particles with a narrow gaussian 
distribution peaked about p = 0), but later consider 
the case where the initial current is non-zero. We take 
p av ~ zky/Dot for particles with positive or negative mo- 
menta respectively; hence \p a v\ estimates the width of 
the momentum distributions about the origin. 

With the approximations < cos2pb >~ 1 and < 
sin2p6 >~ 2pb it is easy to show: 

C ± (2) ~ -K 2 J 2 {K)[cosA{l±b) + 2p av (t)bsmA(l±b)} 

(0.2) 

The ± correspond to the correlations between the two 
possible sequences of two kicks in the map, and we must 
average these to obtain C(2) = 1/2[C(2)_ + C(2)+]. 
However, since b is usually a very small parameter, here 
we can take: 

C(2) ~ -K 2 J 2 (K)[cosA + 2p av (t)bsin A] (0.3) 

or considering separately the energy diffusion for positive 
and negative momenta: 

C(2,±) ~ -K 2 J 2 (K)[cosA±2 v r DolbsmA} (0.4) 

We see immediately from this expression that only the 
terms in sin A differentiate between positive and negative 
momenta and hence allow a ratchet effect. This implies 
that for A = nit, where n = 0, 1, 2.. is an integer, there is 
no transport, while for A = (2n+l)7r/2, there is maximal 
transport. 
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FIG. 1: Figure compares the behavior of the classical current 
as a function of kick strength K for A ~ it/2 , b — 0.03 and 
A ~ 7r. It shows that the latter case is not transporting (ie 
< P >~ 0), but that for A — n /2 the current direction follows 
the oscillations of a Bessel function —J2{K). 

In Fig.l we calculate numerically the classical current, 
as a function of K, for an ensemble of 300,000 trajectories 
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FIG. 2: Comparison between average energy of particles with 
positive momenta E(+) and those moving with negative mo- 
menta E(-) for two different K and b. The figure shows that 
the correction to th e lin ear energy growth is well described at 
short times by Eq. 3.5; at longer times the ratchet correction 
vanishes and E(±) grow as ~ Dot. 



initially sharply peaked about p = (ie a gaussian dis- 
tribution of width a — I). For A ~ tt there is essentially 
no current, while for A ~ 7r/2, the current is apprecia- 
ble and current reversals follow exactly the oscillations of 
— J 2 (K), as expected from the formula. 

In [0 , the origin of the classical current was also found 
to be due to the fact that for p av > we have a (7(2) cor- 
rection of the opposite sign to the correction for p av < 0, 
though in the case of the double- well ratchet the analyti- 
cal form of the correction is much more complicated than 
Eq.0 

In Fig. 2 we test Eq. |0.5| against a numerical calcula- 
tion of the energies as a function of time for K = 3 and 
K = 10 and A = ir/2. For an ensemble of trajectories 
initially peaked about p = 0, we expect that the average 
energy,i? + , of particles with positive momenta, and the 
average energy, _E_ , for those with negative momenta, 
are given by: 

2E± ~ D t ± \bK 2 J 2 (K) VA^ 3/2 (0.5) 

in the regime where pb is small. For large pb, < 
sin2p6 >~< cos2pb >~ hence (7(2) = 0; as shown 
in Q , beyond a time scale t r that it takes the momenta 
to increase to |j3 d^| ~ 1 both the positive p and nega- 
tive p part of the ensemble of trajectories diffuse at the 
same rate 2E =< p 2 >= Dot, where for this system 
D ~ K 2 /2(l - 2J l {K) 2 ..). Hence t r ~ l/(Db 2 ). 

The graph shows that for short times, the equation 
gives a good estimate to the correction to the energy, (< 
p 2 > —Dot). It also shows that after t ~ t r the correction 
vanishes, and < p 2 > —Dot becomes a constant. 

It is very interesting to consider also the case where 
at t = we already have a drift current; in other words 
we start our trajectories with non-zero p = po- It is easy 
to show that the corresponding correction to the energy 



diffusion,C(2,pO)Pati)i now takes the form: 

~ -K 2 J 2 (K) {cos A(cos2p b + 2p± v {t)bsm2p b) 
+ sin A(sin 2p b + 2p± v (t)b cos 2p b)} 

Pa V (t) now represents the average momenta relative to 
Po at time t for particles with momenta greater or less 
than po respecively. We see that the average energy ab- 
sorption rate of the cloud of atoms depends on po as well 
the the effective widths of the distribution p^ v ~ ±^/Dot 
about po- In other words, initially, if all particles are 
at p — po , we have a local diffusion rate correction 
— K 2 J 2 {K) cos{A — 2p b); but as the cloud spreads in 
phase space one has to consider the average over the 
width of the cloud. 

When pb is large, the cos 2pb and sin 2pb terms average 
to zero, hence (7(2, po) — ¥ asymptotically for t > t r re- 
gardless of the value of po- After this, the average energy 
growth for the particles with positive and negative mo- 
menta (relative to p ) is equal to D t. Regardless of po, 
the timescale for these correlations to be 'averaged-out' 
is given by the ratchet time of ||, t r ~ (Dob 2 )^ 1 , the 
time required for the distribution to broaden to a sub- 
stantial width in pb. In H it was found that although 
the classical map does not have the 2-7T periodicity of 
the standard map, it has a new long-ranged periodic- 
ity in p, with periodicity ~ 2n/b. From the formula, 
we see that the periodicity of the correction is half of 
this: (7(2, po = 0) = (7(2, p = 7r/6). We note that since 
the (7(2) correction vanishes when the distribution size 
becomes of order one cell in p, the long range periodic- 
ity in this respect is not significant to the magnitude of 
current accumulated: by the time the expanding distri- 
bution samples the periodicity boundary, < p > is almost 
constant: we have tested this numerically. 

We now consider how one might exploit Eq j0.6| to con- 
struct a cold atom 'turnstile'. For simplicity we take the 
case A = ir/2. Consider the case of a packet of cold 
atoms drifting through the lattice with a constant drift, 
which may be positive or negative, pob — we see that 
for short times (where the equation is valid) its energy 
grows linearly, so < p 2 >~ Dt ~ [D - (±K 2 J 2 (K))]t ~ 

^[1 - 2{J 1 {K)) 2 - (±2J 2 {K))]t. In other words for pos- 
itive drift, the energy absorption is reduced relative to 
the basic diffusion term, while for negative drift it is ac- 
celerated. For K ~ 2 we can cancel the diffusion term 
almost completely so for pob = 7r/4 the packet absorbs 
very little energy. 

Fig. 3a shows the effect on two classical ensembles drift- 
ing through the ratchet with speeds of po = ±27 (for 
an effective turnstile it is not essential to have exactly 
po = ±7r/(4&)) and b = 0.03, K = 1.7. Here D ~ 0.8 so 
for positive p, D ~ 0, while for negative p D ~ 1.6. The 
figure shows the initial distributions as well as the average 
momentum distribution after ~ 125 kicks. We see that 
the effect of the pulsed lattice on the two components 
is drastically different. While the cloud moving right is 
only slightly perturbed, the cloud moving left has been 
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FIG. 3: Shows the effect of a pulsed lattice on a system with 
finite initial current. Figure 3a demonstrates the 'turnstile' 
effect on particles moving through the lattice with initial mo- 
menta po = ±27 ~ 7r/(46). The initial momentum distribu- 
tions are shown, as well as the effect of ~ 125 kicks on the 
final distributions. We see that while particles moving right 
have absorbed relatively little energy, the particles moving left 
have substantially increased average kinetic energy. In Fig 3b 
we show that Eq. ).£ accurately describes the energy absorp- 
tion at short times. The upper figure shows energy growth for 
cos A = 7T and pob — ±7r/4. Energy absorption is asymmetric 
about po and is greater than Dot if po,p a v are of different sign, 
but smaller otherwise. The lower figure corresponds closely 
to Fig3a. Energy absorption rates are symmetric about po 
but are ~ 15 times faster for negative momenta relative to 
positive momenta. 



heated to much larger average kinetic energies. Since for 
this case C(2,po,p^ v ) now has no dependence on p av for 
short times, the particles absorb energy symmetrically 
about po = ±7r/ (46). 

The corresponding quantum case was investigated in 
for the equivalent quantum turnstile to show similar 
behaviour we require only that the quantum break-time 
t* should be of the same order or longer than the duration 
of the classical turnstile (ie ~ 100 kicks for the above 
example, though we can adjust this by varying b). For t < 
t* , the quantum behaviour follows closely the classical 
behaviour. For t > t* , the quantum wavepacket localizes 
and absorbs no more energy, thus 'freezing-in' the ratchet 
asymmetry. For the turnstile, if the lattice region is of 
finite extent, the limiting time could be the time the atom 
spends within the optical lattice if this is less than t* . 



In Fig. 3b we further verify Eq.0.6 by plotting the ac- 
tual energy growth for the case A — tt (upper graph) 
and A = ir/2 (lower graph). For the upper graph, 
only the cos A terms contribute; since we take po — 
±7r/4, we can see from both the formula (C(2,po,p^ v ) ~ 
—2bp^ v sin2pofr) and the numerics that for po,p av of the 
same sign, trajectories absorb energy slower that D$t 
while for po,p av of different sign the converse is true. The 
lower figure corresponds closely to the 'turnstile' shown 
in Fig. 3a since we have sin A = n/2. The figure shows 
that for short times, the energy of particles initially with 
Pa = ±7r/4 grows linearly. However the negative momen- 
tum particles absorb energy ~ 15 times faster than the 
particles with a positive drift. 

In conclusion, we have shown the chaotic ratchet effect 
found in || is generic in character and applies to another 
class of ratchets. We show that by means of a simple an- 
alytical formula the effect can be maximised and use this 
to propose a mechanism for manipulating cold atoms. 
In particular we show that we can have a pulsed optical 
lattice that after ~ 100 kicks has imparted considerable 
energy to particles moving in one direction, while parti- 
cles moving in the opposite direction absorb little energy. 
Since the mechanism is based on fully chaotic dynamics 
rather than a regime with stable islands, no position is 
special so no preparation of the initial state is required. 
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